The matter power spectrum in f(R) gravity 
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Modified gravity can be considered as an alternative to dark energy. In a generalized theory of 
gravity, the universe may accelerate while containing only baryonic and dark matter. We study, 
in particular, the evolution of matter fluctuations in f(R) models within the Palatini approach, 
and find that the resulting matter power spectrum is sensitive to a nonlinear dependence on the 
curvature scalar in the gravitational action. The constraints that arise from comparison to the form 
of the observed matter power spectrum tighten the previous constraints derived from background 
expansion by several orders of magnitude. Models in the allowed parameter space are practically 
indistinguishable from general relativity with a cosmological constant when the backround expansion 
is considered. 



I. INTRODUCTION 

One approach to the current dark energy problem in 
cosmology is to consider the observed acceleration of 
the universe as a consequence of gravitational dynam- 
ics that deviate from general relativity at cosmologi- 
cal scales. This can be achieved by modifying the lin- 
ear proportionality to the curvature (Ricci) scalar in 
the Einstein-Hilbert action to a more general depen- 
dence on curvature invariants, seeQ for a recent intro- 
duction and0 for other applications of extended grav- 
ities in cosmology. Here we will consider models in 
which the gravitational Lagrangian is a nonlinear func- 
tion f(R) of the curvature scalar B, and the deviations 
from general relativity become important at small cur- 
vature. Such cases have been shown to lead to effective 
dark energy 0, Furthermore, we adopt the Pala- 

tini formulation, in which an f(R) gravity is decribed 
by second order field equations]]!, H |j| E3 • The expan- 
sion of the universe has been solved in these models and 
recently also succes fully contrasted with various cosmo- 
logical data sets0,E|, see also^3|. Assuming that the 
correction to the Einstein gravity can be parameterized 
in the form ~ Br, the data excludes the inverse curva- 
ture model (with (3 = —1), but shows a slight preference 
for a non-null correction, although the evidence is not 
compelling 12]. 

In this paper we will investigate further the viability of 
the Palatini approach to nonlinear gravity by extending 
the analysis to the constraints arisin g fr om cosmologi- 
cal perturbations. In a previous work|l4| we were con- 
cerned that the modified matter couplings these mod- 
els feature would, while causing the effective large nega- 
tive background pressure, affect the evolution of inho- 
mogeneities in such ways that the primordial fluctua- 
tions might not form into similar large scale structure 
as observed at the present. We solve cosmology to lin- 
ear order and match the calculated shape of the matter 
power spectrum with the measurements of Sloan Digital 



Sky Survey (SDSS)[l5[. The resulting constraints indeed 
reduce the allowed parameter space into a tiny region 
around the simplest correction to f(B) = R, which is the 
cosmological constant, corresponding to (3 = and thus 
f(B) — B — 2A. For a recent discussion of matter per- 
turbations within the metric formalism of f(B) gravity, 
see [lfj. 

We begin with a brief review of Palatini approach to 
f(B) gravity and its cosmology in section II. The scheme 
for deriving the large scale structure constraints and the 
results are presented in section III. We conclude in section 
IV . 



II. PALATINI APPROACH TO f(R) GRAVITY 
A. Field equations 

We will consider the class of gravity theories repre- 
sented by the action 



S = 



^/[i?(^,f^)]+£ m ( 5 ^,$,. 



(1) 

where n = 8nG. The matter action coupling to gravity 
depends only on some matter fields <&,... and the met- 
ric g^v Our notation emphasizes that the Ricci scalar 
defining the gravitational sector of the action depends on 
the two independent fields, the metric and the connection 
rjg 7 , see also Explicitly, we may write 



B = g^R 



where 



But/ 



(2) 



r™\T ai/ . (3) 



The field equations which follow from extremization of 
the action, Eq.Q), with respect to metric variations, are 
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where the energy momentum tensor is as usually 

2 8{yf^C m ) 



rp{m) — 



(4) 



(5) 
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and we have denned for notational convenience 



F = df/dR. 



(6) 



In general relativity F 
Eq-O gives 



1. Contraction of the field 



RF — 2/ = kT, 



(7) 



an algebraic relation between the trace of the energy mo- 
mentum tensor and the scalar curvature R. 

By varying the action with respect to it is found 
that the connection is compatible with the conformal 
metric 



9fj,v — Fg^ v . 
It follows that the Ricci tensor in Eq.© is 



Rl 



2F 2 



F 



2F 



(8) 



(9) 



R^v being the corresponding tensor associated with the 
metric g^ v . Using this we can write the field Eqs.|@J) ex- 
plicitly in terms of the metric g^ v . The covariant deriva- 
tive V we refer to is associated with the Levi-Civita con- 
nection of this metric. Note also that we have assumed 
that the connection is symmetric, whereas in a more gen- 
eral framework of metric affine gravity ^3 torsion is not 
necessarily absent. 



B. Background cosmology 

We assume a spatially flat FRW universe as the back- 
ground. There the line element is 



ds 2 = -dt 2 + a 2 (t)S i: jdx l dx j , 



(10) 



and matter can be described as a perfect fluid — 
diag(p,p,p,p). Since we will consider late cosmology 
where p — p m = p c + Pb (jn for matter, c for cold dark 
matter, b for baryons), we assume hereafter that p = to 
simplify our equations. The modified Fricdmann equa- 
tion is then 




1 

6F 



(K P + f). 



(11) 



An overdot means derivative with respect to the time 
coordinate t and H is the Hubble parameter H = a/ a. 
By using the conservation of matter and proceeding 
as in^^i we nn d that the Hubble parameter may be ex- 
pressed solely in terms of R, 



H 2 



1 



"if -RF 



6F 



3F' RF- 
2F RF' 



(12) 



For a known function f(R), one can be solve R at a given 
a from Eq. Q and thus get the expansion rate Eq. I|12l) by 
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FIG. 1: The 68, 90 and 99 % confidence contours arising 
from fitting the CMB shift parameter. 



just algebraic means. In what follows, we will also need 
the derivatives H and H. Expressions for these can be 
derived similarly 1 . 

In order to investigate these models quantitatively, we 
must specify the form the gravitational Lagrangian. We 
will use the same parameterization as|l2j. 



f(R) =R-aR f3 1 



(13) 



where a is positive and has dimensions of H 2 ~ 2/3 . The 
(dimensionless) exponent [3 is less than unity, lest the 
correction to the Einstcin-Hilbcrt action would interfere 
with the early cosmology. Given (3 and the amount of 
matter in the present universe, fi m = Kp m /(3H 2 ), deter- 
mines the scale a. Not all combinations of a and are 
consistent with a flat matter dominated universe today, 
and therefore we rather use f2 m and fi as our pair of pa- 
rameters. We will refer to the limit (3 = as the ACDM 
case. 

To demonstrate that models defined by the Lagrangian 
(|13fl can produce a plausible expansion history, one can 
compare the predicted evolution of the Hubble param- 
eter, Ea. l|12fl to cosmological data. As an example we 
compute the CMBR (cosmic microwave background ra- 
diation) shift parameter 19J given by 



dz' 



(14) 



where z,iec is the redshift at decoupling. For these pa- 
rameters we use best-fit values found from the CMBR 



1 However, it is unnecessary to report the resulting (rather 
lengthy) formulae here, since when resorting to numerical anal- 
ysis one may equally well evaluate numerically the derivatives 
of Eq. I12I . the algebraic formulae providing just a means of a 
consistency check. 
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analysis of the WMAP team 20], K = 1.716 ± 0.062 and 
z dec = 1088^2- We plot the constraints resulting from fit- 
ting the CMBR shift parameter in Fig. Projecting to 
the (a,/3) plane would reproduce Fig. 1 of 12], where are 
presented also further constraints derived from multiple 
data sets. 



III. CONSTRAINTS FROM THE LARGE 
SCALE STRUCTURE 

The foundations of cosmological perturbation theory in 
the Palatini approach to modified gravity have been laid 
in|l4|. There was also derived an evolution equation for 
the comoving energy density perturbation 5 m in a matter 
dominated universe. Using the time variable x = log(a), 
the equation may be written as 

^ = A{x ^ + B{x) + JL cll)K , (15) 

where the dimensionless functions A, B and C are given 
by 

FH 2 (2FH + F)A = -2 (h + 2H 2 ^j F 2 

+2FHF + (2FH + FH - 2FH 2 } F, 

FH 4 (2FH + F)B = 2 (H + 2HH^j H 2 F 2 
+2HHFF - (2HHF - HHF - 2H 2 F - 2H 2 HF] F, 



C = ^. (16) 

3F(2FH + F) 

We have also checked that our solutions of the Eq. I|15l) are 
consistent with the explicit field equations listed injl4j]. 

We get the initial conditions for the matter pertuba- 
tion, 

x PcSc + PbSb 

o m = ; , (17) 

Pc + Pb 

by evolving the standard Eistein-Boltzmann equations up 
until zi — 200, corresponding to Xi — — log(201). It is 
only at smaller redshifts that the corrections to general 
relativity begin to have effect. Note also that treating 
the energy density of the universe as consisting of single 
dust-like fluid (made of dark matter and baryons) is well 
justified. At Zi = 200 radiation is subdominant and may 
be neglected, and the imprint it has left on perturbations 
in the earlier universe is carried to the initial conditions 
(namely, 5 m (k,Xi) and dS m (k, Xi)/dx). The Compton 
scattering with photons and pressure of the baryons are 
completely negligible at Zi and therefore baryons obey 
the same equations as dark matter. Hence we can, at the 
linear order, consider just the total fluid. 
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FIG. 2: The matter power spectra for f3 = -0.00001 (solid), 
P = -0.00005 (dash-dotted), (3 = 0.00001 (dotted) and (3 = 
0.00005 (dashed). In all the cases Q, m = 0.3 and ns = 1. (We 
do not use the data points at the three smallest scales, and 
instead of the error bars plotted here we use the exact window 
functions.) 



The effective sound speed j2l| is given by Ea. lfTT))) . and 
is nonvanishing except at the ACDM limit. The curva- 
ture corrections induce effective pressure fluctuations in 
matter, leading to the gradient term in Ea. (|15fl . While 
the other deviations from the perturbation evolution in 
the ACDM scenario may be small when |/3| is, the gradi- 
ent can still be large at small enough scales. As one could 
expect |l4j. inhomogeneities are significantly affected by 
the additional matter couplings, the most sensitive ef- 
fect being the response of modified gravity to spatial 
variations in the distribution of matter. The situation 
bears some resemblance to dark energy models with a 
coupled dark sector [2lL I22I l2flj . In fact, in the confor- 
mally equivalent Einstein frame, where the line-element 
(|10fl is given by the metric (|SJ), one finds a scalar field 
with non- minimal coupling to matter [24|. 

We calculate the matter power spectrum, 

P(k)^(2nk)- 3 S 2 n (k,x = 0), (18) 

where S is considered in the Fourier space. To fit to the 
shape of the observed power spectrum, we use the best- 
fit normalization for each model. We assume that the 
spectral index lies between ns £ [0.8, 1.2] and marginal- 
ize over these values. We make a \ 2 nt to the SDSS 
data 2 0]- Conservatively, we use only measurements at 
scales k < 0.2h Mpc -1 , since at smaller scales there 
are nonlinear effects in the measured P(k). We set the 



2 We use the window functions provided at 
http://space.mit.edu/home/tegmark/sdss.html 
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FIG. 3: The 68, 90 and 99 % confidence contours arising from 
fitting to the SDSS data. The upper panel is for ns = 1 ±0.2 
and the lower panel corresponds to scale-invariant spectrum 
of primordial fluctuations. 



Hubble parameter to h = 0.72 and fix the baryon den- 
sity to fib = 0.044 (since especially the latter is well 
determined 20], and now these parameters have only a 
slight effect on the initial conditions). 

The resulting power spectra is shown for a few choices 
of the parameter /3 in Fig. |3 For these small \(3\, the 
gradient term begins to affect perturbations at just about 
k ~ 0.1 Mpc -1 . For much smaller \/3\, the gradient could 
be neglected at non-linear scales. For much larger \(3\, 
there would be appreciable deviation from ACDM cos- 
mology also at very large scales (due to the ^-independent 
coefficients A and B in the evolution equation Ea. ((T3|l % ). 
and the effect of the gradient would be large enough to 
render these models completely incompatible with SDSS 
data. Thus we can understand why fitting the models to 
the data provides the tight confidence limits plotted in 



Fig. |3 For j3 > 0, red tilt for the primordial spectrum 
is preferred, since the power at small scales is enchanced 
due to modified gravity. Correspondingly, when (3 < 
the best fit is achieved with blue tilt, ns > 1. The prior 
we have assumed, ns € [0.8,1.2], is rather loose, since 
the WMAP experiment is able to rule out large devia- 
tions from scale invariance|20j. at least when standard 
assumptions hold in the pre-recombination universe. For 
this reason we plot also, in the lower panel of Fig. |3J 
the somewhat tighter constraints which follow from fix- 
ing the spectral index to r>g = 1. Marginalizing over 
O m , we find that the favored values of the exponent (3 
are smaller than ~ 3 • 10~ 5 in magnitude. Clearly, inclu- 
sion of the smallest scale SDSS data points would have 
tightened these constraints even further. 



IV. CONCLUSION 

We calculated the matter power spectrum in the Pala- 
tini formulation of modified gravity, and found that the 
observational constraints reduce the allowed parameter 
space to a tiny region around the ACDM cosmology. 
We investigated a specific form for the gravitational La- 
grangian, Eq. (|13|l . but similar conclusions would prob- 
ably hold regardless of the parameterization employed. 
As already found in^4|, the problematic features in the 
perturbation evolution, governed by Eq. (|15l) , are generic 
to these models. 

It would seem difficult to find a way to produce the 
observed large scale structure in the presence of the mod- 
ified matter couplings peculiar to these alternative grav- 
ity theories when F is not very nearly a constant. If this 
could be accomplished by invoking exotic properties to 
the matter components (isocurvature initial conditions, 
or may be entropic perturbations to cancel the effects 
of the couplings), it would not be without adding fine 
tuning and ad hoc assumptions to these models. 

Determining consequences of modified gravity at scales 
of the Solar system[25j or even particle physics[26| exper- 
iments present different challenges as, while exact pre- 
dictions may be incalculable, there can be ambiguities 
in how to parameterize and interpret the possible devia- 
tions from Einstein gravity [271123 ]. Furthermore, it seems 
that by simple arguments broad classes of models can al- 
together evade the Solar system constraints pel l29|. In 
contrast, such implications within cosmological pertur- 
bation theory as discussed here are well understood and 
robust. Hence it might be worthwhile to study to which 
extent the results obtained here can be generalized to 
other forms of alternative gravity theories. This requires 
detailed work, but would enable probing the possibili- 
ties of a more fundamental theory of gravitation by the 
present and up-coming high-precision measurements of 
cosmological large-scale structure. 
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